Modeling  and  Simulation  of  MEMS 

Microthrusters 


Deborah  A.  Levin, 

Associate  Professor, 
Department  of  Aerospace  Engineering 
Penn  State  University 

June  21-25,  2004 

Istanbul,  Turkey 

AFOSR/EOARD  Symposium  on  Energy 
Conversion  Fundamentals 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

21 JUN  2004 

2.  REPORT  TYPE 

N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Modeling  and  Simulation  of  MEMS  Microthrusters 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Department  of  Aerospace  Engineering  Penn  State  University 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADM001793,  International  Symposium  on  Energy  Conversion  Fundamentals  Held  in  Istanbul, 
Turkey  on  21-25  June  2005.,  The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

uu 

18.  NUMBER 
OF  PAGES 

25 

19a.  NAME  OF 
RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


pennState 


MEMS-Based  Propulsion 


•  To  enable  formation  flying  concept  new  propulsion  systems  are 
needed  that  can  deliver  precise  impulse  bits  under  strict  mass,  size 
and  power  limitations. 

•  Various  MEM-based  propulsion  systems  are  considered  for  such 
missions:  cold  and  heated  gas,  bipropellant,  catalytic  and  solid 
decomposition  thrusters. 
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Physics  of  Micro-nozzle  Flows 


•  Due  to  the  reduced  physical  size,  frictional  surface  effects  can  dominate 
the  gas  flow  in  microthrusters. 
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Comparison  of  boundary  layers  for 
traditional  axisymmetric  vs  MEMS  flows. 


Development  of  3D  boundary  layer 
in  a  cold  gas  flow. 


•  In  addition,  wall  temperature,  heat  transfer,  and  heat  fluxes  are  major 
controlling  factors  for  microthruster  performance,  yet  it  is  often  an  unknown 
in  the  system  design. 


The  DSMC  Method 


•A  simulation  tool  for  modeling  chemically  reacting  flows  in 

rarefied! transitional  environments. 

• Kn  =  AAref,  X  =  mean  free  path,  lref  =  reference  length 
— Continuum:  Navier-Stokes,  Kn  <  ~  0.001 
— Transitional:  DSMC  ,  Kn  >  ~  0.01 

•Developed  by  G.  Bird  to  obtain  a  solution  of  the  Boltzmann 
equation, 

—  follow  the  motion  of  many  virtual  molecules  on  a  grid,  for  a 
series  of  time  steps, 

—  calculate  particle  collisions  using  Monte  Carlo  techniques. 

—  model  gas-surface  interactions, 

—  using  conservation  relationships,  obtain  changes  in  internal  energies 
and  velocities  of  the  components. 

•Flows  around  spacecraft  and  micro-propulsion  devices  can  be 
modeled  by  DSMC  because  they  have  similar  Kn. 


Outline  of  T alk 


•  Cold  gas  thrusters  -  effects  of 

-  geometry:  3-D  ,  2-D,  axisymmetric, 

-  gas  surface  model, 

-  performance. 

•  Higher  temperature  gas  thrusters  -  effects  of, 

-  flowfield  dependence  on  T0,  p0 

-  gas-  surface  interaction, 

-  stagnation  conditions  and  internal  energy  on  performance. 

•  High-temperature  gas  thrusters  with  variable  Twam 

-  coupled  DSMC/FEM  method, 

-  variable  material  cooling  conditions, 

-  time-dependent  calculations  of  thrust  and  material 
characteristics. 


Cold  Gas  Thrusters  -  Geometry 


z 


►3D  flat  nozzle 

-XY  plane  expansion  angle  a  =15° 
-throat  width  w  =  300  |um 

-throat  height  h  =  300  jum 


-Area  ratio 


Ae/A*  =  10 


•Axisymmetric 

-expansion  angle  a  =15° 
-throat  radius  Rt  =  150  jum 


-Area  ratio 


Ae/A*  =  100 


•2D  nozzle,  neglect  surfaces  in  x-y  plane,  for  z=0,  h 


Cold  Gas  Thruster  -  Conditions 


•Test  gas 

•Stagnation  temperature  and  pressure 
•Critical  temperature  and  pressure 
•Wall  temperature 

•Knudsen  number  (mean  free  path/char  length) 
•Reynolds  number  (mom/viscous  ratio) 


N2 

Tc  =  300  K,  Pc=  10  kPa 
Tt  =  250  K,  P,=  5.2  kPa 
Tw = 300  K 
Kn  =  5x1  O'3 
Re  =  200 
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Cold  Gas  -  Effect  of  Geometry 

Top  -  3-D,  Bottom  -2-D 
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•Density,  kg/m3 
•Greater  flow 
expansion  in  2-D 
case. 


•T,  K 

•More  heat  transfer 
to  the  wall  in  the  3- 
D  case. 
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Cold  Gas  -  Effect  of  Geometry 

3-D  vs.  Axisymmetric 
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DSMC,  3D 
DSMC,  AXISYM 


Solns  are  the 
same  since 
mass  flow 
rates  equal. 


0.0075 


Solns  are  the  same  since  mass  flow 


•3-D  vs  axisymmetric  shows  that  the  degree  of  wall  surface  area  is  important. 
•For  supersonic  nozzle  flow  into  a  vacuum,  Ux  should  be  a  maximum  at  the  exit. 
•Instead,  extremum  Ux  is  located  upstream  due  to  subsonic  region  at  the  walls. 
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Kinetic,  DSMC  Gas-Surface  Wall  Models 

in  Micro-nozzle  Flows 


Specular  reflection  Diffuse  reflection 


incident  gas 


reflected  gas 
molecule  „ 


•Two  types  of  gas-surface 
interactions 

-Specular  -  incident  and 
reflected  tangential 
momentum  are  the  same, 

-Diffuse  -  all  reflected 
directions  are  equally 
probable,  wall  acts  like 
an  emitting  source  of 
particles  at  Twall, 


•Interaction  is  specified  by  with  accommodation  coefficients  for 

oNormal  momentum,  tangential  momentum,  and  energy  (heat  transfer) 
oValues  between  0,1  and  closer  to  1  for  spacecraft  materials. 

•Continuum  methods  are  partially  corrected  with  wall  slip-jump  boundary  conditions. 
•Accommodation  ~  1  generates  large,  viscous  boundary  layers. 
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Cold  Gas  Thrusters  -  Effect  of  Gas- 

surface  model 


■W -  DSMC,  alpha-d  =  0. 


■e -  DSMC,  alpha-d  =0.3 
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Translational  T  profiles  along  the  nozzle  axis  for 
different  ad  in  an  axisymmetric  micronozzle. 


•Definition  of  tangential 
momentum  accommodation 
coefficient, 

ad=(P,i  -  J\, )/  Pxi 
•Experiments  show  that  ad  = 
0.8  for  silicon. 

•3-D  calculations  showed 
only  1%  difference  for  ad  = 
0.8  and  1. 
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Cold  Gas  Thrusters  -  Performance 
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Case 

Thrust  (mN) 

Up,  (sec) 

AS  NS 

1.07 

65.62 

AS  DSMC 

1.03 

65.5 

2D  NS 

1.17 

69.45 

2D  DSMC 

1.10 

68.74 

3D  DSMC 

0.93 

56.61 

X,  m 

•Comparison  of  3-D  calculations  with  data  show  good  agreement. 

•2-D  assumption  is  poor  and  over  predicts  thrust  levels  of  geometry. 

•Wall  effects  in  3-D  case  reduce  thrust  (20%)  and  specific  impulse  cf  to  2-D. 
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Higher-Temperature  Thruster  Flows 

Full  wall  accommodation,  Ret  ~  200,  Tw  =  300  K,  3-D  nozzle 


Normalized  density 


Ux  along  centerline 


•Flows  at  both  temperatures  are  dominated  by  surface  interactions,  but  structure 
is  different. 

•Isp=  56.6  and  61.s  for  To=300, 1000  K,  respectively. 

•Increase  in  Isp  « than  for  comparable  axisymmetric  case  due  to  larger  surface- 
area-volume  ratio. 
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Comparison  of  Flows  for  Different  Gas  - 
Surface  Models  and  T0 


1.  “Specular”  =  ideally  smooth,  no  momentum  or  energy  transfer  (TOP), 

2.  “Diffuse,  adiabatic”  =  av  tangential  momentum  of  reflected  molecules  =  0,  no  energy 
transfer  with  wall  (MIDDLE). 

3.  “Diffuse,  Tw=300”  =  both  momentum  and  energy  transfer  occurs  with  wall 
(BOTTOM). 
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Effect  of  T0  and  P0  on  Performance 

Axisymmetric  nozzle,  Diffuse  wall,  Tw=300 


f/f,  KP> sec  To~2,300  K 


•For  high  T0,  a  shorter  nozzle  would  give  better  performance  at  lower  Re. 
•Need  to  optimize  geometry  of  high  temperature  nozzles. 

•For  lower  pressures,  the  peak  value  of  Isp  is  closer  to  the  nozzle  throat. 
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Effect  of  Internal  Energy*  on  Performance 


Axisymmetric  nozzle,  hydrogen 
air  mixture  66.1%  N2, 32.4% 
H20,  and  1.5%  H2 
High  T  polyatomic  gas,  VT 
becomes  important. 

Represent  VT  by  Zv  relaxation 
numbers  in  Larsen-Borgnakke 
model: 

1.  104  <  ZV(T)  <  106 

2.  Constant  Zv  =  1  and  100. 
Faster  VT  relaxation  increases 
nozzle  performance. 


*For  cold  gas  thrusters,  TR  relaxation  is  the  dominant  internal  energy  transfer  mechanism. 
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Comparison  of  Calculated  and  Measured 
Efficiencies  for  Axisymmetric  Nozzles 


Efficiency,  % 


•  Experimental  data  did  not 
specify  Tw. 

•  Good  agreement  between 
measurements  and  modeling 
for  Tw=500  K  and  complete 
accommodation. 

•  Ln  sensitive  to  wall 
conditions. 

•  Need  a  predictive  capability  to 
determine  both  material  and  gas 
properties  for  micropropulsion 
devices . 
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Geometry  and  Flow  Conditions 
Schematic  of  NASA-Glen  Microthruster 


30  deg  converging  part 
15  deg  diverging  part 
Exit  to  throat  area  ratio  of  5 
Throat:  300pm  x  600  pm 
Nitrogen  flow: 

-  P0  =  0.1  and  0.5  atm, 

-  To=2000  K, 

.  _  -  Re  =  35  and  175, 

1.7  mm 

respectively. 


10  m 


8.8  mm 
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Coupled  Thermal  and  Gas  Dynamic 
Computational  Approach 


0.0)5 
0.0)4 
0.003 
0.002 
0.001 
0 

Computational  domain  for  2D  calculations 

•  Based  on  the  solution  of  the  heat  transfer  problem  using  finite  element 
method  coupled  to  the  DSMC  gas  flow  solution. 

•  Coupling  between  material  thermal  response  and  flow  by  using  DSMC 
heat  duxes  as  the  boundary  conditions  for  the  heat  conduction  problem. 

•  The  wall  temperature  calculated  in  the  heat  transfer  simulations  is,  in 
turn,  applied  as  a  boundary  condition  for  DSMC  calculations. 
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Material  thermal  response:  3D 

Time  variation  of  the  temperature  fields  for  0.1  atm 


Thermally  Insulated 


Active  Cooling 
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Micronozzle  performance 
h=600  pm,  po=0.1  atm,  insulated  vs  cooling 


•Final  thrust  value,  F,  are  15%  and  6%  lower  than  the  initial 
ones  for  Case  1  and  2,  respectively. 

•The  mass  flow  degradation  is  as  much  as  55%  in  Case  1. 

•Coeff.  of  mass  discharge,  Cd,  decreases  more  rapidly  than  F. 


Where  does  the  energy  go? 


•2-D,  po=0.5  atm, 
cooling,  (AIAA  03-0673). 

•Cooling  allows  heat 
transfer  losses. 

•As  material  heats  up, 
less  gas  energy  is 
transferred  to  the  wall. 

•As  material  heats  up, 
reduction  in  thrust  will 
slow  down. 


Ratio  of  gas  energy  to  value  at  inlet 
along  the  nozzle  axis. 


Conclusions 
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Material  thermal  response  is  similar  to  2D,  but  gas  flow  structure  is 
significantly  different  due  to  the  side-wall  boundary  layer. 

For  thermal  cooling,  the  steady-state  material  temperature  is  450  K  (same  as  in 
2D). 

Higher  Re  flow  results  in  larger  surface  heat  fluxes  and,  thus,  shorter 
operational  times.  However,  heat  fluxes  do  not  vary  proportionally  with  Re. 

Cooling  applied  to  the  outer  surface  can  sustain  the  material  temperature 
below  melting.  Cooling  results  in  improved  thrust  and  mass  discharge 
performance. 

The  large  temporal  variation  of  the  thrust,  and  especially  mass  discharge 
coefficient  means  that  the  coupling  between  the  gas  and  material  must  be  taken 
into  account  in  micropropulsion  design. 
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